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FOREWORD 


This  report  describes  an  in-house  assessment  by  the  author  for  the 
High  Speed  Aero  Performance  Branch,  Aeromechanics  Division,  Air  Force 
Wright  Aeronautical  Laboratories  (AFWAL/FIMG) .  The  work  was  performed 
under  Project  2404,  "Aeromechanics,"  Task  24040730,  "Aerodynamic 
Prediction  Methods." 


The  evaluation  reported  herein  used  methodology  developed  by  the 
Aerophysics  Research  Corporation,  authored  by  F.  A.  Woodward  and  by  the 
McDonnell  Douglas  Corporation,  authored  by  A.  E.  Gentry,  D.  N.  Smith,  and 
W.  R.  Oliver.  The  experimental  data  necessary  to  evaluate  these  method¬ 
ologies  were  extracted  from  technical  reports  published  by  the  Air  Force 
Armament  Laboratory,  authored  by  C.  B.  Butler,  E.  S.  Sears,  and  S.  G. 
Pallas;  NASA,  authored  by  E.  J.  Landrom;  and  AEDC,  authored  by  D.  C.  Baker 
and  D.  E.  Reichenau.  The  evaluation  reported  herein  was  conducted 
between  October  1979  and  July  1980. 
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SECTION  I 
INTRODUCTION 

An  engineer  is  often  faced  with  j  difficult  task  in  selecting  the 
analysis  method  that  will  provide  good  and  consistent  results.  One  of  the 
things  that  complicates  the  selection  is  model  geometry  complexities,  such 
as  irregular  body  shapes  and  control  surfaces  that  lead  to  complex  flow 
fields  beyond  the  analysis  capabilities  of  the  more  simplified  analytical 
procedures.  A  relative  comparison  between  analytical  methods  should  be 
made  with  simple  geometries  so  that  under  these  relaxed  conditions  com¬ 
parisons  between  various  methods  might  lead  to  some  insights  into  their 
relative  capabilities  and  restrictions.  Thus  for  the  present  study, 
slender  ogive/cylindrical  body  combinations  without  control  surfaces, 
nozzles,  etc.,  were  selected  to  study  the  relative  merits  of  analysis 
methods  that  the  Flight  Dynamics  Laboratory  often  uses  to  analyze  more 
complex  configurations.  Low  supersonic  Mach  numbers  were  selected  for 
this  study  as  well  as  a  subsonic  case  because  experience  has  shown  that 
these  two  cases  are  the  most  difficult  to  predict.  Selected  for 
evaluation  in  this  study  were  the  Unified  Subsonic-Supersonic  Aerodynamic 
Computer  Program  (USSAERO)  of  Reference  1  and  the  Supersonic-Hypersonic 
Arbitrary  Body  Program  (HABP)  of  Reference  2. 

The  objective  of  the  present  study  is  to  take  a  detailed  look  at  the 
capabilities  of  the  above  programs  for  predicting  overall  aerodynamic 
characteristics.  Such  things  as  detailed  pressure  coefficient  predictions 
around  the  body  are  compared  with  the  experimental  data  to  determine  where 
the  analysis  would  benefit  from  improvement.  Configuration  geometry 
variables  such  as  nose  bluntness  and  fineness  ratio  are  varied  in  the 
supersonic  case  to  study  their  impact  on  the  analysis.  Comparisons  are  made 
at  a  subsonic  Mach  number  of  0.4  and  supersonic  Mach  numbers  beginning  with  a 
low  value  of  1.50  and  progressively  increasing  to  4.6.  The  comparisons  in 
the  present  study  have  been  restricted  to  the  merits  of  the  inviscid 
methods  because  inclusion  of  the  viscous  analysis  techniques  would  have 
complicated  the  evaluation  and  perhaps  led  to  confusion  in  interpreting 
the  results  contained  herein. 


AFWAL-TR-81 -3035 


SECTION  II 
MODEL  DESCRIPTIONS 

1 .  N2B1  (FROM  MX  TESTS) 

A  large  number  of  missile  configurations  were  used  in  the  present 
investigation  for  the  data/theory  comparisons.  Complete  descriptions  of 
these  models  and  their  instrumentation  can  be  found  in  References  3,  4, 
and  5.  A  brief  description  is  repeated  here  for  completeness.  For  the 
subsonic  comparisons,  the  configuration  was  the  N2B1  configuration 
(Reference  3).  The  model  had  an  ogive  nose  with  a  fineness  ratio,  LfJ/D, 
of  2.545  fineness  ratio  and  an  aft  body  of  7.754  fineness  ratio  for  a 
combined  fineness  ratio  of  10.3.  The  bluntness  ratio  R^/R^,  was  .03. 
Figure  1  shows  a  sketch  of  the  model  geometry  and  pressure  orifice 
locations . 

2.  CAC  AND  BNC  (FROM  NASA  TESTS) 

For  the  supersonic  pressure  prediction  comparisons,  two  models  from 
Reference  4  were  selected.  They  are  the  circular-arc-cylinder  (designated 
CAC  in  the  present  study)  model  and  the  Lil  unt-nose-cyl inder  (designated 
BNC)  model,  with  a  bluntness  ratio  of  0.17.  The  geometrical  description 
of  these  two  configurations  are  as  follows: 

Circular-Arc-Cylinder  (CAC) 

r/L  =  '  (1.3125)2  +  (X/L) (0.9  -  X/L)  -  1.3125  (0  <  X/L  <  0.45) 

r/L  =  0.075  (0.45  <  X/L  <  1 .0) 

Blunt-Nose-Cyl inder  (BNC) 

r/L  =  0.1118  (X/L)1/2  (0  <  X/L  <  0.45) 

r/L  =  0.075  (0.45  <  X/L  <  1.0) 

A  sketch  of  these  two  models  showing  model  orifice  locations  is  seen  in 
Figure  2. 


AFWAL-TR-81 -3035 


3.  AFATL  MODELS 

For  the  force/moment  comparisons  at  supersonic  speeds  a  large 
selection  of  configurations  were  available  from  tests  conducted  by  the 
Air  Force  Armament  Test  Laboratory  (AFATL)  as  reported  in  Reference  5. 
The  models  selected  for  the  present  comparisons  consisted  of  20  nose 
geometries  with  fineness  ratios  varying  from  1  to  4  and  one  a  ft  body 
with  a  fineness  ratio  of  10.  Bluntness  ratios  (R^/Rg)  varied  from  sharp 
0.0  to  1.0.  The  geometrical  description  and  designations  of  the  models 
are  presented  in  Figure  3. 
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SECTION  III 
METHODOLOGY  OVERVIEW 


1.  UNIFIED  SUBSONIC-SUPERSONIC  AERODYNAMICS  PROGRAM  (USSAERO) 

A  complete  description  of  the  theory  and  computer  program  may  be 
obtained  from  Reference  1 .  Presented  here  is  a  brief  summary  of  the 
general  theory  used  for  body  alone  applications.  Axis  and  coefficient 
nomenclature  is  presented  in  Figure  4. 

The  configuration  surface  is  subdivided  into  a  large  number  of 
panels,  each  of  which  contains  an  aerodynamic  singularity  distribution. 

A  constant  source  distribution  per  unit  area  is  used  on  each  body  panel. 
The  position  of  the  source  within  each  panel  relative  to  a  coordinate 
system  lying  in  the  plane  of  the  panel  is  denoted  by  (c.,  n,  o).  The 
radius  connecting  an  arbitrarily  located  source  within  the  quadrilateral 
with  a  point  in  the  flow  field  P(X,  Y,  Z)  is  given  by 

h  =  C(x  -  T  )2  +  (y  -  r,}2  +  z2]1/2  (1) 


for  incompressible  flow  applications.  Details  for  the  compressible  flow 
case  can  be  found  in  Reference  1  and  will  not  be  summarized  here.  The 
general  approach,  however,  is  to  apply  Gothert's  rule  to  the  incompressible 
velocity  components  described  in  this  section.  The  incremental  perturbation 
potential  arising  from  an  increment  of  area  (df.d-i )  located  at  the  source 
point  U,  n)  is 


Thus 


d^dn 


(2) 


$  = 


Jll 

4tt 


// 

A 


cMn 


C(X  -  cr  +  (y  -  o /  + 


7p 


(3) 


The  general  approach  of  the  method  is  to  calculate  the  velocity  components 
at  the  centroid  of  each  panel  by  partial  differentiation  of  the  potential 
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defined  in  Equation  3  with  respect  to  x,  y,  z.  The  pressure 
coefficient  at  the  centroid  of  the  jth  panel  is  then  calculated  as 
fol lows 


where 

q2  =  (V  +  cosi)2  +  V  2  +  (Vz  +  sin,.)2 
J  j  j  j 

The  total  force  and  moment  coefficients  are  then  obtained  by  summing  the 
products  of  the  pressure  coefficient  over  the  entire  body  geometry. 

Details  of  the  mathematics  required  to  obtain  the  velocity  components  on 
each  of  the  panels  are  shown  for  the  incompressible  case  in  Section  I  of 
the  Appendix. 

2 .  SUPERSONIC-HYPERSONIC  ARBITRARY  BODY  PROGRAM  ( HABP) 

The  Mark  IV  Supersonic-Hypersonic  Arbitrary  Body  Program  (HABP)  is  a 
computer  program  that  is  capable  of  calculating  the  aerodynamic  character¬ 
istics  of  arbitrary  complex  three-dimensional  shapes.  A  complete 
description  cf  the  program  and  its  various  pressure  coefficient  options  is 
provided  in  Reference  2.  Although  the  program  is  most  accurate  at  hyper¬ 
sonic  speeds,  its  capabilities  have  been  extended  down  to  supersonic  speeds 
by  use  of  slender  body  theory  and  empirical  techniques  that  are  described 
in  this  section.  As  with  USSAERO  the  body  is  described  by  a  large  number 
of  panels.  Element  area,  centroid,  and  a  unit  normal  are  then  calculated 
from  the  four  points  describing  the  quadri 1 ateral  surface.  Before  the 
program  calculates  the  pressure  on  each  surface  element,  it  checks  to  see 
if  the  element  is  facing  the  flow  (in  an  impact  region).  This  is  done 
by  checking  if  the  sign  of  the  angle  n  between  the  freestream  velocity 
and  the  outward  directed  unit  normal  vector  to  the  surface  is  positive. 

If  the  sign  is  negative,  the  selected  expansion  method  is  used.  The  unit 
outward  directed  normal  on  each  of  the  surface  elements  is 

-*  -» 

n  =  nxi  +  ii  j  +  n2k  (5) 
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where  n  ,  n  ,  and  n  are  the  direction  cosines  of  the  surface  element 
x  y  z 

normal  necessary  to  obtain  the  components  of  the  force  on  each  quadrilateral. 
The  mathematical  details  for  obtaining  these  direction  cosines  is  outlined 
in  Reference  2. 


A  common  parameter  necessary  to  all  of  the  inviscid  methods  in  HABP 
is  the  angle  6  that  a  tangent  to  the  element  makes  with  the  freestream 
velocity.  Once  the  unit  outward  directed  normal  expressed  in  Equation  5 
is  obtained  on  a  particular  quadrilateral,  the  tangent  angle  is  determined 
simply  as  follows : 


where 


(6) 


There  are  15  inviscid  compression  pressure  methods  and  nine  expansion 
pressure  methods  available  as  an  option  to  HABP  users.  The  selection  of 
the  proper  method  in  a  given  application  depends  upon  the  vehicle  component 
and  flight  condition  and  is  selected  on  the  basis  of  knowledge  and 
experience  in  the  use  of  each  method. 

For  the  present  study,  three  combinations  of  compression/expansion 
inviscid  pressure  methods  have  been  selected  for  comparison  with  data  on 
the  lower  end  of  the  Mach  number  range  (Mach  1.5  to  4.63).  They  are  as 
follows : 

Compression  Methods 

1  -  Newtonian 

5  -  Approximate  Cone  (Slender  Body  plus  Bammitt-Murthy) 

14  -  Modified  Dahl  em/Buck 

Expansion  Methods 

1  -  Newtonian  (C  =  0. ) 

P 

3  -  Prandtl  -Meyer 

5  -  Van  Dyke  Unified 
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In  the  comparisons  with  the  data,  a  combination  of  a  compression  method 
and  an  expansion  method  is  designated  by  a  two-number  sequence.  The 
first  number  represents  the  impact  method  and  the  second  number  the 
expansion  method. 


a.  Newtonian  Compression  Method 


The  Newtonian  theory  on  the  impact  surfaces  gives  the  pressure 
coefficient  as  a  direct  function  of  the  tangent  angle  '  (Equation  6)  as 
fol lows 

Cp  =  Ksin26  (7) 
where  y 


K 


max 


Y  +  3 

Y  +  1 


[1 


2 

Y  +  3 


] 


b.  Approximate  Cone  Compression  Method 

This  method  is  used  in  HABP  to  improve  the  low  Mach  number 
pressure  coefficient  prediction  capabilities.  Three  approximate  cone 
techniques  are  used.  Specifically,  second-order  slender  body  theory 
(Reference  6)  is  used  in  the  low  supersonic  range  and  the  technique 
described  in  Reference  7  is  used  in  the  higher  supersonic  range.  In 
addition,  an  empirical  tangent  cone  method  described  in  Reference  2  is 
used  for  large  incidences  where  the  flow  becomes  detached. 

The  second-order  slender  body  solution  of  Reference  6  is  used 
whenever  (Mmsin6)  <  a  where  a  =  0.2  for  freestream  Mach  numbers  less  than 
3.0  and  a  =  0.3  for  larger  freestream  Mach  numbers.  Then  from  Reference  6 
the  pressure  coefficient  is  a  function  of  5  and  M^  and  is  determined  as 
fol lows : 

Cp  *  tan26{ (2c  -  1.0)  +  tan26 [3(8C )2  -  (5M2  -  1 )c  +  A])  (8) 

cone(2nd  order) 

where  <5  is  determined  as  in  Equation  6  and 

f  .  2 
^  in  Btane 
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For  the  higher  Mach  number  the  Hammi tt-Murthy  solution  of  Reference  7  is 
used  whenever  (M^sinS)  >  0.325.  In  this  method  the  shock  cone  angle, 
os,  and  surface  Mach  number,  M  ,  is  first  determined.  Results  from 
Reference  8  are  then  used  to  obtain  the  pressure  ratio  and  local  Mach 
number  behind  the  shock.  Isentropic  relations  between  the  shock  and 
body  are  then  used  to  obtain  the  pressure  ratio  on  the  surface,  Pc/P  . 
The  pressure  coefficient  is  then  found  as  follows: 


l'orie  (Hammi  tt-Murthy) 


(Pr/P  )  -  1.0 

HM  _ 

1/2yM2 


(9) 


The  mathematical  details  of  the  Hammi tt-Murthy  solution  are  presented 
in  Section  2  of  the  Appendix.  If  the  flow  is  detached  a  different 
procedure  is  used  to  obtain  o  .  Details  are  also  found  in  Section  2 
of  the  Appendix.  If  (M  sins)  is  between  0.3  and  0.325,  a  transition 
solution  between  the  second-order  results  and  the  Hammitt-Murthy 
solution  is  determined.  The  details  are  found  in  Section  2  of  the 
Appendix. 


c.  Modified  Dahl  em/Buck  Compression  Method 

For  high  Mach  number  (>  20)  this  method  uses  an  empirical 
relationship  that  approximates  tangent-cone  pressures  at  low  impact 
angles  and  approaches  Newtonian  values  at  high  impact  angles.  The 
relationship  is 


where 


and 


Cp  =  Ksi n 


M  >  20 


K  = 


1  .0 


+  1.0 


(sin46)0,75 
*  2.0  for  6  >  22.5° 


for  6  <  22. 5C 


(10) 
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For  small  values  of  6  the  upper  limit  of  the  bracketed  term  is  not 
allowed  to  exceed  5.0.  The  method  is  extended  to  lower  Mach  numbers 
by  assuming  that  the  empirical  pressure  coefficient  of  Equation  10  varies 
with  Mach  number  in  the  same  manner  as  the  pressure  coefficient  for  a 
right  circular  cone. 

This  gives 


C 


PMDB 


=  Cr 


(Eq  10) 


Pcone(Mi  •  20) 
Pcone(Mu  =  20) 


(11) 


The  term  in  brackets  in  Equation  11  was  determined  empirically  in 
Reference  9  to  be 


where 


and 


CP 

cone(M  <  20) 

- =  —  ,  -S  in  degrees 

bp  r  n 

'cone^  =  20) 


a  =  (6.0  -  0.3M  )  +  sin 

cjo  ' 


0 . 588ti 


) 


n  =  1.15  +  0.5  sin 


(  tn(M J  -  0.916^ 
\  3.29  / 


(12) 


d.  Newtonian  Expansion  Method 

This  method  substitutes  a  lower  limit  on  the  expansion  pressure 
equal  to  the  freestream  value.  Thus,  for  6  <  0  the  pressure  coefficient, 
Cp,  is  set  equal  to  zero. 

e.  Prandtl -Meyer  Expansion  Method 

This  method  utilizes  the  relationships  of  Reference  8  to  expand 
the  freestream  Mach  number  to  a  surface  Mach  number,  M  .  With  Mc  thus 
obtained  the  familiar  isentropic  relationships  of  Reference  8  are  then 
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used  to  calculate  the  ratio  of  the  surface  pressure  to  total  pressure. 
The  pressure  coefficient  is  then  calculated  as  follows: 


vX> 


P  /PT 
c  T 

P../P  t 


1.0 


(13) 


r/(y  -  1) 


-  1.0 


\  M 

2  ' 


(14) 


where 


rd  m2 

7  -  (i  +  V  w~) 
rlM2 

O  +  2  Mc) 

Details  of  the  expansion  process  are  found  in  Section  3  of  the  Appendix. 

f.  Van  Dyke  Unified  Expansion  Method 

This  method  for  expansion  flow  is  derived  in  Reference  10. 

The  results  are 


- 

'P  v  2 

VAN  DYKE  yH 


(1.0 


2y/(y-D 

— ) 


1.0 


for  H  >  2/(y-1)  (15) 


where 


H  =  (i/m2  -  TTo)< 


and 


Cp  =  *  4  for  H  <  2/U-l) 

vAN  DYKE  M 
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SECTION  IV 

COMPARISONS  OF  THE  ANALYSIS  WITH  EXPERIMENT 

It  should  be  noted  that,  in  all  of  the  comparisons  to  follow,  base 
pressure  has  been  subtracted  out  and  skin  friction  has  not  been  included 
in  the  USSAERO  comparisons  or  the  integrated  pressure  HABP  results. 

1 .  USSAERO  MX  DATA 

Selected  portions  of  the  data  base  from  Reference  3  were  utilized 
to  compare  the  USSAERO  capability  to  predict  the  detailed  pressure 
coefficients  and  forces  on  typical  ogive  nose  shapes.  The  comparisons 
of  the  analysis  and  experiment  are  described  here  for  the  nose  geometry 
shown  in  Figure  1  and  for  Mach  0.4  freestream  flow. 

Figure  5  shows  the  comparison  of  USSAERO  with  measured  pressure 
coefficients  at  20  degrees  angle  of  attack  and  the  total  integrated 
normal  coefficients  for  an  angle  of  attack  range  from  0  to  20  degrees. 

For  model  stations  between  X/D  =  0.46  and  1.84  the  theory  compares  very 
well  with  the  measured  results  for  pressure  coefficients.  On  the  remaining 
20  percent  of  the  nose,  however,  USSAERO  predicted  negative  pressure 
coefficients  on  the  windward  meridians  and  positive  pressure  coefficients 
on  the  leeside.  This  results  in  the  prediction  of  large  negative  loads 
on  the  aft  end  of  the  nose  with  a  resulting  under-prediction  of  the 
total  normal  force  coefficient  as  evidenced  on  the  right  in  Figure  5. 
Included  for  comparison  purposes  is  the  simple  expression  from  Reference  11 
that  relates  the  normal  force  to  angle  of  attack  and  the  crossflow  drag 
coefficient.  As  seen  from  the  figure,  the  inclusion  of  the  crossflow 
drag  component  improves  the  correlations.  Analytical  studies  using  a 
source  distribution  approach  similar  to  Woodward's  method  (References  12, 
13)  experienced  the  same  problem  on  short  axisymmetric  bodies.  These 
studies  attributed  the  source  of  the  problem  to  the  lack  of  proper 
modeling  of  the  body  wake.  The  correction  used  in  these  studies  was  to 
model  the  base  streamlines  as  part  of  the  body  description.  This  was 
done  by  extending  the  body  and  then  closing  it  with  a  wake  of  the  same 
shape  as  the  nose. 
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In  the  present  study,  the  effect  of  extending  the  body  downstream 
without  closure  was  made  to  determine  its  effect  on  the  Woodward  analysis. 
Pressure  coefficient  data  was  available  from  Reference  3  for  cylindrical 
body  extensions  to  the  nose  geometry  of  Figure  1.  The  nose  had  a  fineness 
ratio  of  =  2.55.  Cylindrical  body  extensions  of  0.25  diameters 

and  2.25  diameters  were  added  to  the  N2  nose  producing  two  additional 
configurations  with  total  fineness  ratios  of  2.C  and  4.8,  respectively. 

An  additional  comparison  was  made  with  the  N2  with  a  7.75  diameter  long 
afterbody  added  to  produce  a  configuration  with  a  total  fineness  ratio  of 
10.3.  The  result  seen  on  the  left  side  of  Figure  6  shows  the  effect  of 
the  small  aft  body  extension  (fineness  ratio  =  2.8)  and  the  larger 
extension  (fineness  ratio  =  4.8)  on  the  Woodward  predictions  at  the 
nose/body  juncture  (X/D  =  2.3)  and  on  the  extended  cylindrical  body. 

As  noted  in  the  figure  the  addition  of  the  2.25  body  extension  corrected 
the  negative  load  prediction  at  the  base  of  the  N2  nose  (denoted  as  the 
nose/body  juncture  in  this  figure).  In  implementing  this  procedure 
another  problem  in  the  Woodward  analysis  became  apparent.  This  can  be 
seen  on  the  left  side  of  Figure  6  for  the  pressure  coefficient  comparisons 
on  the  body.  Note  that  while  the  angle  of  attack  is  20  degrees  the 
Woodward  analysis  predicts  a  nearly  symmetrical  pressure  distribution  on 
the  extended  nose/body  combination  of  L/D  =  4.8;  thus  resulting  in  zero 
normal  force  on  the  cylinder.  This  problem  can  most  likely  be  attributed 
to  the  absence  of  a  solution  to  the  2-D  Laplace  equation  in  the  Y,  Z 
crossflow  plane  in  the  Woodward  analysis.  Without  this  addition  one 
might  not  expect  good  agreement  with  data  at  angles  of  attack  greater 
than  about  four  degrees. 

This  is  indeed  the  case  as  evidenced  on  the  right  of  Figure  6 
where  lack  of  agreement  for  angles  of  attack  above  four  degrees  become 
significant  for  longer  bodies  as  the  comparisons  with  the  10.3  fineness 
ratio  configuration  show. 

Thus,  the  present  study  shows  that  (1)  for  a  small  fineness  ratio 
nose,  extending  the  nose  downstream  one  or  two  diameters  will  correct 
the  negative  loads  on  the  aft  end,  however,  (2)  when  attempting  to 
predict  the  loads  on  nose/body  combinations,  the  absence  of  a  crossflow 
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solution  in  the  subsonic  Woodward  analysis  at  low  Mach  numbers  may 
result  in  significant  errors  in  predicting  the  normal  force  coefficient 
at  angles  of  attack  above  four  degrees  (especially  for  large  fineness 
ratio  nose/body  combinations).  No  additional  analysis  was  made  to 
determine  if  closing  the  base  with  an  extension  of  the  same  shape  as  the 
nose,  such  as  done  in  References  1?  and  13,  would  improve  the  correlations. 

2.  HABP-NASA  PRESSURES 

The  remaining  studies  of  this  report  compare  supersonic  predictions 
using  various  combinations  of  the  HABP  compression  and  expansion  methods 
outlined  in  Section  III  and  in  the  Appendix.  The  comparisons  are  made  for 
the  two  nose/body  combinations  shown  in  Figure  2  and  for  the  large 
variety  of  nose  shapes  shown  in  Figure  3.  Results  are  presented  over  a 
number  range  between  1.63  to  4.63  for  the  nose/body  combinations  and 
consist  of  comparing  the  predictions  with  the  detailed  pressure  coefficient 
data  from  Reference  4.  For  the  various  nose  shapes  seen  in  Figure  3,  the 
comparisons  consist  of  total  force  and  moment  comparisons.  This  data  was 
obtained  from  Reference  5  at  Mach  numbers  from  1.5  to  4.0. 


The  particular  HABP  compression/expansion  prediction  combination 
used  in  the  comparisons  are  identified  in  all  of  the  remaining  figures  as 
a  two-digit  number  (c,  e)  where  c  denotes  the  compression  method  and  e 
the  expansion  method.  The  following  numbers  have  been  selected  to 
identify  these  selections. 


Compression 


c  =  1 
c  =  5 
c  =  14 


Newtonian 
Approximate  Cone 
Modified  Dahlem/Buck 


Expansion 

e  =  1 
e  =  3 
e  =  5 


Newtonian  (Cp  =  0) 
Prandtl -Meyer 
Van  Dyke  Unified 


Comparisons  for  the  pressure  coefficient  predictions  are  presented  in 
Figures  7  through  19.  The  solid  and  dashed  lines  show  the  theoretical 
predictions.  The  circles  represent  the  experimental  data  from  Reference  4. 
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To  be  consistent  with  the  data  format  in  Reference  4  the  nomencl atore  of 
the  circumferential  location,  ,  differs  from  that  presented  in  Figures  5 
and  6  of  this  report.  The  nomenclature  of  in  Figures  7  through  17  is 
as  shown  in  Figure  2.  o  =  0°  is  the  centerline  on  the  upper  surface  and 
0  =  180°  is  the  lower  surface  centerline. 

Each  figure  is  arranged  to  show  relative  effects  of  Mach  number  and 
nose  bluntness  at  a  selected  angle  of  attack  for  a  particular  X/L  location 
corresponding  to  the  instrumentation  stations  in  Figure  2.  For  example, 
in  Figure  7a  comparisons  between  the  blunt  and  sharp  model  are  shown  on 
the  left  of  the  figure  for  Mach  1.6  while  the  same  comparisons  are  made 
on  the  right  for  Mach  4.63.  All  of  these  comparisons  are  made  for  X/L 
values  from  0.075  to  0.225.  Identical  comparisons  are  shown  for  the 
remaining  axial  stations  in  Figures  7b  through  7d.  All  comparisons  in 
Figure  7  are  for  alpha  equal  to  -4  degrees.  Figure  8  shows  the  distributed 
loads  obtained  by  integrating  the  experimental  and  predicted  pressure 
distributions  from  Figure  7  along  the  body.  Figures  9  through  18  repeat 
this  pattern  for  the  remaining  angles  of  attack  through  alpha  -  16  degrees. 
Finally,  the  total  integrated  loads  at  each  of  the  angles  of  attack  are 
shown  in  Figure  19. 

A  careful  study  of  the  pressure  comparisons  show  that  methods 
(5,3)  and  (14,5)  are  generally  more  accurate  on  the  nose  (X/L  between 
0.125  and  0.425)  while  method  (1,1)  is  in  better  agreement  on  the 
cylindrical  aft  body.  On  the  cylinder  the  expansion  prediction  Cp=0 
is  in  general  agreement  with  the  data.  The  overall  correlations  improve 
with  increasing  Mach  number.  The  resulting  total  force  and  moment 
characteristics  as  shown  in  Figure  19  reflect  that  method  (1,1)  provides 
better  agreement  with  and  because  they  are  more  influenced  by 
pressure  prediction  accuracy  on  the  cylinder.  The  remaining  two  methods 
have  more  accuracy  for  the  predictions  since  the  axial  force  is  more 
influenced  by  the  nose  pressure  prediction  accuracy.  Figure  19  shows  that 
methods  (5,3)  and  (14,5)  were  more  accurate  in  predicting  for  the 
sharp  model  than  for  the  blunt  model.  However,  since  these  comparisons 
are  only  for  one  bluntness  ratio  (R^/Rg  =  0.17)  and  one  fineness  ratio 
(L^/D  =  3.2),  an  expanded  data  base  comparison  featuring  variations  in 
these  variables  will  be  presented  in  the  remaining  paragraphs. 
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3.  HABP-AFATL  FORCE  AND  MOMENTS 

The  data  base  from  Reference  5  was  selected  for  the  geometrical 
comparisons  because  of  the  large  selection  of  nose  bluntness  and  fineness 
ratios  available  as  seen  from  Figure  3.  The  unit  Reynolds  number  for 
this  data  was  R  / ft  =  1.8  x  10^  for  all  Mach  numbers.  The  boundary  layer 
was  assumed  to  be  turbulent  and  the  reference  enthalpy  skin  friction 
option  was  selected  in  the  remaining  HABP  predictions.  The  predicted 
skin  friction  coefficient  variation  was  about  the  same  for  all  of  the 
configurations  analyzed  and  varied  nominally  from  0.04  at  Mach  1.5  to 
0.03  at  Mach  4.0.  Since  the  skin  friction  coefficient  was  nearly  the 
same  for  each  of  the  nose  shapes  at  a  particular  Mach  number,  the 
relative  comparisons  of  predicted  are  sensitive  to  the  inviscid 
method  selected.  The  comparisons  in  the  following  figures  are  made  from 
the  experimental  data  with  the  base  pressure  drag  subtracted.  Since  most 
of  the  ogive  nose  surface  experiences  a  compression  in  the  angle  of 
attack  range  from  0  to  16  degrees,  the  remainder  of  the  study  focused  on  the 
HABP  compression  methods.  The  Newtonian  expansion  method  (Cp  =  0)  was 
arbitrarily  selected  for  each  combination  of  compression  methods  and 
remained  fixed  for  all  of  the  geometrical  comparisons. 

The  remaining  force  and  moment  comparisons  are  shown  in  Figures  20 
to  26.  Each  figure  is  arranged  to  show  the  effect  of  nose  bluntness 
and  Mach  number  for  a  particular  fineness  ratio.  For  example,  Figure  20a 
shows  comparison  between  the  sharp  nose  and  a  nose  with  bluntness  ratio 
of  1.0.  The  nose  fineness  ratio  is  also  1.0.  Comparisons  on  the  left 
are  for  Mach  1.50  while  those  on  the  right  of  the  figure  are  for  Mach  2.0. 
Figure  20b  continues  these  comparisons  for  Mach  3.0  and  Mach  4.0.  The 
remaining  figures  follow  the  same  pattern  of  presentation  for  successively 
increasing  fineness  ratios.  For  the  higher  fineness  ratio  nose  shapes, 
a  larger  selection  of  nose  bluntness  data  was  available  for  comparison: 
generally  ranging  between  R^/Rg  =  0.25  to  0.75. 

A  study  of  the  comparisons  shows  that  each  of  the  compression  methods 
has  satisfactory  agreement  with  the  data  for  bluntness  ratios  less  than 
U.4,  even  down  to  Mach  numbers  as  low  as  1.5.  For  R^/Rg  =  1.0  and  0.75 
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each  of  the  compression  methods  generally  over-predicted  the  axial  force 
at  all  Mach  numbers.  For  smaller  bluntness  ratios,  this  was  not  always 
true,  especially  with  the  Newtonian  compression  method.  With  one 
exception,  Mach  1.5,  each  of  the  three  methods  under-predicted  axial 
force  for  the  sharp  nose  models  at  all  Mach  numbers.  By  comparing  the 
sharp  model  results  for  progressively  higher  fineness,  one  can  note  a 
slight  improvement  in  predictions  as  fineness  ratio  increases.  The 
bluntness  ratio  is  the  dominant  geometrical  variable  affecting  the 
analysis  capability,  and  the  agreement  improves  for  all  of  the  methods 
as  the  Mach  number  increases.  The  Newtonian  compression  method  (C  =  1) 
was  the  most  accurate  in  predicting  the  normal  force  and  pitching  moment 
coefficients  over  the  entire  Mach  number  range;  this  was  the  same 
result  obtained  in  the  nose/body  configuration  studies  (Figures  7  through 
19). 
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SECTION  V 
CONCLUSIONS 

An  evaluation  has  been  made  on  two  aerodynamic  analysis  programs 
that  one  might  use  to  predict  the  aerodynamic  characteristics  of  simple 
slender  bodies  without  complexities  such  as  control  surfaces,  nozzles, 
etc.  Specifically,  in  the  subsonic  range  USSAERO  (Reference  1)  which 
subdivides  the  body  into  a  large  number  of  panels  with  a  constant 
source  distribution  on  each  panel,  was  se’ected  for  evaluation.  In 
the  low  supersonic  regime  (M  =  1.5  to  4.63)  three  impact  methods  and 
three  shadow  methods,  available  in  the  Supersonic-Hypersonic  Arbitrary 
Body  Program,  referred  to  as  HABP,  were  studied.  From  this  study  a 
number  of  conclusions  were  reached  and  are  presented  as  follows: 

The  USSAERO  program  gave  satisfactory  subsonic  predictions  on  the 
ogive  nose  geometries  for  angles  of  attack  less  than  eight  degrees. 

For  higher  angles  of  attack,  two  problems  were  encountered.  The  first 
problem  was  a  large  negative  load  predicted  on  the  aft  end  of  the  nose. 
This  was  a  result  of  negative  pressure  coefficients  on  the  windward  rays 
and  positive  pressure  coefficients  on  the  leeward  rays  being  predicted. 
The  problem  was  corrected  by  extending  the  nose  2.25  diameters  downstream 
with  a  cylindrical  afterbody.  In  implementing  this  procedure,  a  second 
problem  in  the  USSAERO  analysis  became  evident.  The  second  problem  was 
the  prediction  of  zero  normal  loads  on  the  cylindrical  extensions  even 
at  angles  of  attack  up  to  20  degrees.  This  was  a  result  of  predicted 
symmetrical  pressure  distributions  on  the  cylindrical  extensions.  It  is 
felt  that  this  is  due  to  the  20  crossflow  plane  solution  not  being 
present  in  the  subsonic  USSAERO  analysis. 

Three  compression  and  three  expansion  HABP  prediction  methods  were 
evaluated  from  Mach  1.5  to  Mach  4.63  for  a  large  variety  of  ogive  nose 
shapes.  The  primary  variables  were  nose  bluntness  and  fineness  ratio. 

In  addition  comparisons  were  made  for  a  nose/body  combination  of  overall 
fineness  ratio  of  6.67  with  both  a  sharp  and  blunt  nose.  The  results  for 
the  nose/body  combination  showed  that  approximate  cone  and  modified 
Dahlem/Buck  compression  predictions  methods  in  conjunction  with 
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Prandtl -Meyer  and  Van  Dyke  Unified  expansion  methods  were  more  accurate 
than  the  Newtonian  predictions  for  due  to  the  larger  contribution  of 
the  nose  on  the  predictions.  On  the  other  hand,  the  Newtonian  methods 
were  superior  for  the  and  predictions  which  are  influenced  more  by 
the  cylindrical  afterbody  than  the  nose. 

Comparisons  for  a  large  variety  of  ogive  nose  shapes,  with  nose 
bluntness  and  fineness  ratio  being  the  primary  variables,  showed  that 
nose  bluntness  was  the  significant  geometrical  variable  that  determined 
prediction  accuracy.  For  bluntness  ratios  above  R^/Rg  =  0.4  all  of  the 
methods  over  predicted  C^.  The  sharp  nose  exhibited  slight  sensitivity 
to  nose  fineness  ratio  in  predicted  C^.  The  normal  force  coefficient, 
C^,  and  the  pitching  moment  coefficient,  C^,  predictions  were  relatively 
insensitive  to  these  geometrical  variables. 
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APPENDIX 

DETAILED  FORMULATION  OF  METHODS 

Some  of  the  intermediate  procedures  necessary  in  obtaining  the 
pressure  coefficients  in  Section  III  are  included  here  for  completeness 

1.  USSAERO/ DISTRIBUTED  SOURCES 

As  seen  in  Equation  4  the  pressure  coefficient  acting  at  the 
centroid  of  each  panel  is  determined  by  the  components  of  velocity  in 
the  X,  Y,  and  Z  direction.  These  are  obtained  by  the  partial  differen¬ 
tiation  of  the  source  potential  defined  in  Equation  3  with  respect  to 
the  x,  y,  and  z  coordinates  referenced  to  the  leading  of  each  panel  as 
shown  in  the  following  sketch.  The  resulting  velocities  are  then  trans 
ferred  to  the  components  in  the  X,  Y,  Z  coordinate  system. 

To  begin  the  procedure  consider  the  geometry  of  each  panel  as 
defined  by  the  following  sketch 
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The  velocity  components  in  the  x,  y,  z  directions  are  obtained  by 
differentiating  Equation  3.  The  resulting  expressions  are: 


_x  _  jl  rr 

V  ^x  ~  4t  j/ 
a)  A 


(x-rjdc  dn 


[(x-c)2  +  (y-n)2  +  z2]3/2 


(A-l ) 


-1  =  *  =  _L  rr  _ (y-n)dc_dn _ 

V  y  4*  r/  .2  .  .2  A  2,3/2 

°°  A  [(x-f,)  +  (y-n)  +  z  J 


( A-2) 


JL  =  *  .  JL  ;; 

*„  *2  A 


Z  dcdn 


C(x-C)2  +  (y-n)2  +  z2]3/2 


(A-3) 


Complete  details  of  the  derivation  may  be  found  in  Reference  1.  However, 
the  result  for  a  typical  corner  is  included  here  for  completeness.  With 
a  =  tan6  (where  <5  is  found  as  in  Equation  6)  and  iiu  =  the  body  panel  edge 
slope  dy/dx,  the  three  velocity  components  at  each  of  the  four  corner 
points  are  determined  as  follows  in  the  X,  Y,  Z  coordinate  system. 


- : - ?~T7T  K  G  -  H  -  aF],  i  =  1  to  4 

wi  i 


_  -GO  ,  a2)1''2 

4tt 


=  1  to  4 


- : - PHT?  +  a  (miG  -  H)]>  i  =  1  to  4 

4n(l  +  a2)  ' 2  1 


( A-  4 ) 


( A- 5 ) 


( A- 6 ) 
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where 


F  =  tan 


'( 


(z  -  ax) ( 


2  2  2  ^  ^  \ 

l  *r  *  n . ) 

-  z  (ay-m.z)  / 


-x  (y-mi  x)  -  z  (ay-m.z) 


G  = 


1.0 


(1  +  a2  +  m2)1/2 


sinh’^ 


(x  +  m.y  +  az  \ 

[(y-mix)2+(ay-miZ)2+(z-ax)2]1/2/ 


H  =  si 


1nh’1  ((x2  *  77*) 


The  velocity  components  at  each  of  the  four  corner  points  are  found 
by  applying  the  above  formulas  with  the  origin  shifted  to  the  corner  under 
consideration  and  using  the  appropriate  edge  slope.  The  influence  of  the 
panel  is  obtained  by  summing  the  influences  of  the  four  corners  and  applying 
the  resulting  velocity  to  the  centroid  of  the  panel.  Thus,  for  the  jth 
panel 


V 


j 


vY  ) 

4  j 


+  V 


( A- 7 ) 
( A— 8) 
( A- 9 ) 


These  velocity  components  are  then  substituted  into  Equation  4  to  deter¬ 
mine  the  pressure  coefficient  at  the  centroid  of  the  panel. 

2.  HABP/HAMMITT-MURTHY 

The  objective  of  the  Hammitt-Murthy  solution  is  to  determine  the 
surface  pressure  ratio,  Pc/P^,  for  application  in  Equation  9  to  calculate 
the  surface  pressure  coefficient.  The  procedure  used  is  to  determine  the 
shock  cone  angle,  os,  and  surface  Mach  number  Mc>  as  outlined  in 
Reference  7.  The  pressure  ratio  and  local  Mach  number  behind  the  shock 
are  next  obtained  followed  by  an  isentropic  relation  between  the  shock 

and  surface  to  obtain  P  /P  . 
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To  determine  0s  the  basic  equation  to  be  solved  is,  from  Reference  7, 
the  quadratic 

H1  h2  +  (2  sin  26)  h$  +  H£  =  0  (A-10) 

where 

hs  -  (0S  -  6) 

H]  =  2  -  (Y  +  5)  sin2  6 
H2  =  - C( Y - 1 )  sin26  +  2/M2] 

Note  that  once  hg  has  been  calculated  and  6  has  been  found  from  Equation  6, 
C>s  is  then  determined.  Equation  A-10  is  first  checked  to  see  if  the 
roots  are  real  or  imaginary.  If  the  roots  are  imaginary  the  flow  is 
detached  and  the  Tangent  Cone  Empirical  method  of  Reference  2  is  used, 
as  fol iows : 

sino  =  [K  M  sin  6  +  e'KcM°»  Sin6]/M  (A-ll) 

where 

Kc  =  2(y+l )/(y+3) 

If  the  roots  of  Equation  A-10  are  real  then  the  solution  is  taken  as 
the  quadratic  root 

-npx  H.  H~  1/2 

0=6-  ^—2.  1.0  -  [1.0  +  - 5 — - —  ]  (A-12) 

S  H1  sin2  (26) 

If  H^O,  then  the  term  within  the  brackets  in  Equation  A-12  is  expanded 
into  a  power  series  and  the  solution  becomes 

es  “  5  *  2sTnT2sT  -0  "  T  H1H2  0-0  -  ^  (A-13) 
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The  surface  Mach  number  is  determined  as 

flf  cos2qs  (l+2h2) 

'  c  1  +  ^  Mf  (sin2os  -  2h2  cos20$) 

The  pressure  ratio  behind  the  shock  is  determined  next  from  Reference  8 

P  /P  =  4r  M2  sin20  -  *4  ( A-l 5) 

S  00  ^+1  00  S  Y+  1 

Next  the  isentropic  relations  of  Reference  8  are  used  between  the  shock 
and  surface  for  constant  Ty.  First,  the  surface  temperature  ratio  is 
determined. 


1/2 


(A-l 4] 


]c__ 

T-  (1.0  +I§iM2) 

The  ratio  T  / T  is  found  as  follows  (Reference  8) 

°o  s 


(A-l 6) 


where 


(Ps/Poo  +  A) 

Ps  Ps 

jr-  (1-0  +  A  p^-) 


oo 


CO 


(A-l 7) 


Next  the  isentropic  relation 


of  Reference  8  is  used  to  determine  Pc/Ps 


Pc/P 


s 


( A-l  8) 


Finally,  the  surface  pressure  ratio  P  /P^  is  determined  as 


P 

oo 


(A-l 9) 
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~  is  then  substituted  into  Equation  9  to  obtain  the  surface  pressure 

oo 

ratio.  It  was  pointed  out  in  the  main  text  that  the  Hammi tt-Murthy 
solution  is  applied  for  M^sind  >  0.325  and  the  second  order  solution  for 
M  sind  <  0.3.  For  values  of  M  sind  between  these  values  a  transition 

oo  oo 

solution  is  obtained.  This  is  found  by  using  a  Taylor's  series  expansion 
of  the  second  order  solution  in  the  neighborhood  of  M^sind  =  0.3. 

Specifically  let 

Cp(x)  =  Cp  for  a  <  M^sind  <  b  (A-20) 

C0NE(TRANSITI0N) 

where  a  =  0.3  and  b  =  0.325 


and 


Cp(a)  =  Cp  for  M  sind  =  a 

C0NE(2nd-order) 

Cp(b)  =  Cp 

C0NE( HAMMI TT-MURTHY) 

X  =  M  sind,  AX  =  X-a 

00 

h  =  b-a 


(A-21) 

(A-22) 

( A- 23) 
( A-24 ) 


Then  the  transition  solution  is  found  as 

dC  d2C  2 

cp(x)  -  Cp ( a )  +  Ax  +  2  Ax 

dX 

The  first  derivative  is  found  by  di fferentiating  Equation  8  with 
respect  to  tand  and  applying  the  following  chain  rule: 


(A-25) 


dCP  =  dCP  d(tand) 
d(Mtosindy  ~  d(tandl  ’  d^sindT 

^ P  d(tand )  dd 

d(tand )  dd  ’  d(Moosind) 


dCp 

d(tan<5 ) 


2, 
sec  d 

M  COSd 


dCr 


d{tand ) 


1 

M  C0S3d 


(A-26) 
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Upon  differentiating  Equation  8  with  respect  to  tanfi  and  combining  with 
Equation  A-26  the  result  is 

dCp  P  P 

=  tan6{4(rj  -  1.0)  +  tan  6[6'  (2c  -  1):' 

( A- 2  7 ) 

-  ( 5Mf  -  1 ) (4f>1 )  +  4A]}/(McC0S3  ) 

where  the  terms  in  this  equation  are  the  same  as  defined  in  Equation  8. 
The  second  derivative  is  found  from  the  finite  difference  relation 


d2Cr 


dX 


Cp  (b)  -  2Cp  (a)  +  Cp  (a-h) 


where 


Cp  (a-h)  = 


dC 

Cp  U)  -  jf  •  h 


2  dCp 

d^Cp  Cp  (b)  -  Cp- (a)  -  h 

dX^  li 


(A-28) 


(A-29) 


3.  HABP- PRANDTL -MEYER  EXPANSION 

The  general  procedure  for  the  Prandtl -Meyer  expansion  from  free- 
stream  is  as  follows:  first,  the  relations  from  Reference  8  are  used  to 
obtain  the  surface  Mach  number,  Mc-  Then  an  isentropic  relationship  is 
applied  between  the  freestream  and  surface  to  obtain  PC/POT-  This 
expansion  procedure  is  used  whenever  5  from  Equation  6  is  negative. 

The  Prandtl -Meyer  angle  corresponding  to  the  freestream  is 


v»  tan  1  '  tan  'ft7' 

The  Prandtl -Meyer  angle  after  expansion  is 

vc  =  +  I6! 


( A-30) 


( A- 31 ) 
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The  approximate  downstream  surface  Mach  number  is  found  as  follows 


dMco 

Differentiating  Equation  A-30  for  Au  and  substituting  into 

Equation  A-32  gives 


( A- 32 ) 


where 


(to)O  *  Zj1  M*) 

+7^  . 


( A-33) 


Av  =  \j  -  v 
C  ® 


)«l 


Equation  A-30  can  be  used  with  vc  replacing  and  Mc  replacing  Mr  and 
rearranged  to  give 


Mc  - 


1  .0  +(  R 


Ill 


-1 


1/2 


(A -34) 


where  R  =  tan 


Equation  A-34  has  the  form 


vc  +  tan 


"c  -  S(MC) 


or 


and 


Mc  -  g(Mc)  =  F(M^)  =  0 


(A-35) 


1  -  g’(Mc)  -  F'(MC) 


( A- 36 ) 
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Equation  A-35  may  be  solved  by  iterating  for  the  surface  Mach  number  Mc 
using  the  Newtonian  method 


F(MCj) 

F'(*'  ) 

Li 


Mr  -  g(M  ) 
i  ci 


(A-37) 

( A-38) 


g(M  )  -  g1 (M  )M 
ci  i  i 


(A-39) 


With  Equation  A-32  used  as  a  starting  value,  successive  iterations  of 

Equation  A-39  will  yield  the  final  surface  Mach  number  Mc  for  a  given 

value  of  v  .  With  M  determined,  the  pressure  coefficient  is  determined 
c  c 

from  the  isentropic  relationship  given  in  Equation  14. 
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MODEL  GEOMETRY 


1  Geometry  and  Pressure  Orifice  Locations  for  USSALRO 
ari  sons 


riODEL  PRESSURE  ORIFICE  LOCATIONS 
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Figure  2.  Model  Geometries  an  Pressure  Orifice  Locations  for  HABP 
Pressure  Prediction  Comparisons. 
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1  MOMENT  REFERENCE  LOCATION  AT 
THE  NOSE  CYLINDER  JUNCTURE 

2.  REFERENCE  AREA  = 

3  REFERENCE  LENGTH  =  D 
4.  ALL  DIMENSIONS  IN  INCHES 


USSAERO  BODY  GEOMETRY  REPRESENTATION 


Figure  4.  Axis  and  Coefficient  Nomenclature 


Figure  6.  Effect  of  Adding  an  Afterbody  on  the  USSAERO  Comparisons  with 
Experiment. 


Figure  7.  Comparison  of  Pressure  Coefficient  Data  with  the  HABP  Analytical 
Results.  Effects  of  Mach  Number  and  Nose  Bluntness, 

Alpha  =  -4  Degrees. 


Figure  7.  Continued 


Figure  7.  Continued 


Figure  7.  Concluded 
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Figure  8.  Integrated  HABP  Comparisons  with  Experiment  for  Distributed  Local 
Force  and  Moment  Results,  Alpha  =  -4  Degrees. 


Figure  9.  HABP  Pressure  Comparisons  at  Alpha  =  0  Degrees. 
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Figure  9.  Continued 
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Figure  9.  Continued 
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Figure  9.  Concluded 
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Figure  10.  Integrated  HABP  Local  Force  and  Moment  Distribution, 
Alpha  =  0  Degrees. 


Figure  11.  HABP  Pressure  Comparisons  for  Alpha  -  4  Degrees. 


Figure  11.  Continued 


Figure  11.  Continued 


Figure  12.  Integrated  HABP  Local  Force  and  Moment  Distribution, 
Alpha  =  4  Degrees. 
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Continued 


Figure  13.  Concluded 


Figure  ?4.  Integrated  HA BP  Local  Force  and  Moment  Distribution,  Alpha  =  8  Degrees. 


Pressure  Prediction  Comparisons  at  Alpha  =  12  Degrees. 


Figure  15.  Continued 


Figure  ? 5 .  Concluded 


Figure  16.  Integrated  HABP  Local  Force  and  Moment  Distribution, 
Alpha  =  12  Degrees. 


Figure  17.  Continued 


Figure  17.  Continued 


Figure  17.  Concluded 


Figure  18.  Integrated  HABP  Local  Force  and  Moment  Distribution, 
Alpha  =  16  Degrees. 


Total  Force  and  Moment  HABP  Comparisons  with  Experiment  at 
Various  Angles  of  Attack.  Effect  of  Mach  Number  and  Nose 
Bluntness. 
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Figure  20.  Effect  of  Nose  Bluntness  and  Mach  Number  on  Total  Force  and 
Moment  HABP  Comparisons.  Nose  Fineness  Ratio  =  1.0. 
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Figure  21.  HABP  Force  and  Moment  Comparisons  for  the  2.0  Nose  Fineness 
Ratio  Geometry.  Effects  of  Nose  Biuntness  and  Mach  Number. 
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Figure  21.  Continued 


Figure  21 .  Continued 


Figure  21 .  Continued 


Figure  21.  Concluded 


Figure  22.  Continued 


Figure  24.  HABP  Force  and  Moment  Comparisons  for  the  Power  Series  Nose 
Geometry. 
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Figure  24.  Concluded 


Figure  25.  HABP  Force  and  Moment  Comparisons  for  the  3.0  Fineness  Ratio  Nose. 


Figure  25.  Continued 
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Figure  25.  Continued 


Moment  Comparisons  for  the  4.0  Fineness  Ratio  Nose 


Figure  26.  Continued 
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ire  26.  Continued 
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